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t "divide -and- conquer"  [1].  In  this  method,  the  original  problem  is 

divided  into  two  or  more  smaller  problems.  The  subproblems  are  solved 
by  applying  the  method  recursively,  and  the  solutions  to  the  subproblems 
are  combined  to  give  the  solution  to  the  original  problem.  Divide-and- 
conquer  is  especiall.y  efficient  ■V'^ien  the  subproblems  axe  substantially 
smaller  than  the  original  problem. 

I In  [lit]  the  following  theorem  is  proved. 

I 

Theorem  1.  Let  G be  any  n-vertex  planar  graph  with  non-negative 
vertex  costs  summing  to  no  more  than  one.  Then  the  vertices  of  G can 
be  partitioned  into  three  sets  A,  B,  C , such  that  no  edge  joins  a 
i vertex  in  A with  a vertex  in  B , neither  A nor  B has  total  vertex 

cost  exceeding  2/3  , and  C contains  no  more  than  2'\/2'/n  vertices. 
Furthermore  A,  B,  C can  be  found  in  0(n)  time. 

In  the  special  case  of  equal-cost  vertices,  this  theorem  becomes 


Corollary  1.  Let  G be  any  n-vertex  planar  graph.  The  vertices  of  G 
can  be  partitioned  into  three  sets  A , B , C , such  that  no  edge  joins  a 
vertex  in  A with  a vertex  in  B , neither  A nor  B contains  more 
than  2n/3  vertices,  and  C contains  no  more  than  2'\/2Vn  vertices. 

Theorem  1 and  its  corollary  open  the  way  for  efficient  application 
of  divide-and-conquer  to  a variety  of  problems  on  planar  graphs.  In  this 
paper  we  explore  a number  of  such  applications.  Each  section  of  the 
paper  describes  a different  use  of  divide-and-conquer.  The  results  range 
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from  an  efficient  approximation  algorithm  for  finding  maximum  independent 
sets  in  planar  graphs  to  lower  bounds  on  the  complexity  of  planar 
Boolean  circuits.  The  last  section  mentions  two  additional  applications 
whose  description  is  too  lengthy  to  be  included  in  this  paper. 
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2.  Approximation  Algorithms  for  NP-Complete  Problems. 

Divide-and-conquer  in  ccmbination  with  Theorem  1 can  be  used  to 
rapidly  find  good  approximate  solutions  to  certain  NP-complete  problems  on 
planar  graphs.  As  an  example  we  consider  the  maximum  independent  set 
problem,  which  asks  for  a msiximum  number  of  pedrwise  non-adjacent 
vertices  in  a planar  graph. 

Theorem  2.  Let  G be  an  n-vertex  planar  graph  with  non-negative  vertex 
costs  summing  to  no  more  than  one  and  let  0 < £ 1 • Then  there  is  some 

set  C of  0(  Vn/e  ) vertices  whose  removal  leaves  G with  no  connected 
component  of  cost  exceeding  e . Furthermore  the  set  C can  be  found 
in  0(n  log  n)  time. 

iTOof . Aprly  the  following  algorithm  to  G . 

Initialization;  Let  C = 0 . 

General  Step:  Find  some  connected  conponent  K in  G minus  C with 

cost  exceeding  c . Apply  Corollary  1 to  K , producing  a partition 
Aj^  , of  its  vertices.  Let  C = C U • If  one  of  A^^  and  B 

(say  A^  ) has  cost  exceeding  two-thirds  the  cost  of  K , apply 
Theorem  1 to  the  subgraph  of  G induced  by  the  vertex  set  , 
producing  a partition  , Bg  f Cg  of  A^  . Let  C = C U C2  • 

Repeat  the  general  step  until  G minus  C has  no  component  with 
cost  exceeding  e . 

The  effect  of  one  execution  of  the  general  step  is  to  divide  the 
comronent  K into  smaller  components,  each  with  no  more  than  two-thirds 
the  cost  of  K and  each  with  no  more  than  two-thirds  as  many  vercloes 


as  K . Consider  all  components  which  arise  during  the  course  of  the 
algorithm.  Assign  a level  to  each  component  as  follows.  If  the 
component  exists  when  the  algorithm  halts,  the  component  has  level  zero. 
Otherwise  the  level  of  the  component  is  one  greater  than  the  maximum 
level  of  the  components  formed  when  it  is  split  by  the  general  step. 

With  this  definition,  any  two  components  on  the  same  level  are 
vertex-disjoint. 

Each  level  one  component  has  cost  greater  than  e , since  it  is 
eventually  split  by  the  general  step,  it  follows  that,  for  i > 1 , 
each  level  i component  has  cost  at  least  (5/2)^”^€  and  contains 
at  least  (3/2)^  vertices.  Since  the  total  cost  of  G is  at  most  one, 
the  toteil  number  of  components  of  level  i is  at  most  (2/3  )^~^/  e . 

The  total  running  time  of  the  algorithm  is  0(Z/  {jK]  | K is  a component 
split  by  the  general  step])  . Since  a component  of  level  i contains 
at  least  (3/2)^  vertices,  the  maximum  level  k must  satisfy 
(3/2)  < n , or  k < log^^2  n . Since  components  in  each  level  are 

vertex-disjoint,  the  total  rtuming  time  of  the  algorithm  is 
0(n  logj/g  • 

The  total  size  of  the  set  C produced  by  the  algorithm  is  bounded  by 
o(E  I K is  a component  split  by  the  general  step}) 


^1.1063/2  nj  L(2/3)^'VeJ 

< 0l  S maxV  Z vnl  1 T n,  < n 

\ i = l j = l J ' j = l J - 


nj  < n and  n^  ^ 0 


0(Vn/€  ) . □ 


The  following  algorithm  uses  Theorem  2 to  find  an  approximately 
maximum  independent  set  I in  a planar  graph  G = (V, E)  . 


Step  1.  Apply  Theorem  2 to  G with  € = (log  log  n)/n  and  each  vertex 
having  cost  l/n  to  find  a set  of  vertices  C containing 
O(n/Vlog  log  n ) vertices  whose  removal  leaves  no 
connected  component  with  more  than  log  log  n vertices. 

Step  2.  In  each  connected  component  of  G minus  C , find  a maximum 
independent  set  by  checking  every  subset  of  vertices  for 
independence.  Form  I as  a union  of  maximum  independent  sets, 
one  from  each  component. 


Let  I*  be  a maximum  independent  set  of  G . The  restriction  of 
I*  to  one  of  the  connected  components  formed  when  C is  removed 
fron  G can  be  no  larger  than  the  restriction  of  I to  the  same 
component.  Thus  |l*|  - |ll  = o(n/Vlog  log  n ) . Since  G is  planar, 

G is  four-colorable,  and  |l*|  >n/h  . Thus  (|l*l-|l|)/  ll*|  = 

0(l  /vlog  log  n ) , and  the  relative  error  in  the  size  of  I tends  to 
zero  with  increasing  n . 

Step  1 of  the  algorithm  requires  0(n  log  n)  time  by  Theorem  2. 
n. 

Step  2 requires  0(n.  2 ) time  on  a connected  component  of  n^  vertices. 

The  total  time  required  by  Step  2 is  thus 

^ J n n.  n ^ "N 

0 max<  Tj  n.  2 ^ 1 n.  = n and  0 < n.  < log  log  n \ 1 = 

V V = ^ " i = i " " j j 

0^  log  log  n n)2^°®  ^ = 0(n  log  n)  . Hence  the  entire 

algorithm  requires  0(n  log  n)  time. 
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Nonserial  Dynamic  Programming. 


Many  NP-complete  problems,  such  as  the  maximum  independent  set 

problem,  the  graph  coloring  problem,  and  others,  can  be  formulated  as 

nonserial  dynajnic  programming  problems  [ 2,20].  Such  a problem  is 

of  the  following  form:  minimize  the  objective  function  f(x^,  , 

where  f is  given  as  a sum  of  terms  > each  of  which  is  a function 

of  only  a subset  of  the  variables.  We  shall  assume  that  all  variables 

X.  take  on  values  from  the  same  finite  set  S , and  that  the  values 
1 

of  the  terms  given  by  tables.  Associated  with  such  an 

objective  function  f is  an  interaction  graph  G = (V, E)  , containing 

one  vertex  v^  for  each  variable  x^  in  f , and  an  edge  Joining  x^ 

and  X . for  any  two  variables  x.  and  x which  appear  in  a common 
D id 

term  • 

By  trying  all  possible  values  of  the  variables,  a nonserial 

dynamic  programming  problem  can  be  solved  in  2^^^^  time.  We  shall 

show  that  if  the  interaction  graph  of  the  problem  is  planar,  the 

problem  can  be  solved  in  2*^^*'^^  time.  This  means  that  substantial 

savings  are  possible  ^dien  solving  typical  NP-complete  problems  restricted 

to  planar  graphs.  Note  that  if  the  interaction  graph  of  f is  planar, 

no  term  ^ contain  more  than  four  variables,  since  the 

complete  graph  on  five  vertices  is  not  planar. 

Ill  order  to  describe  the  algorithm,  we  need  one  additional  concept. 

m 

The  restriction  of  an  objective  function  f = S to  a set  of 


variables  x^  ,...,x^  is  the  objective  function 

1 J 

f = S {fu  1 fi,  depends  upon  one  or  more  of  x.  , ...,x  ] 

K K Ij 
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Given  an  objective  fiinction  f(x  , ...,x  ) = £ f i,  and  a 

^ ” k=  1 ^ 

subset  S of  the  variables  x^, ...,x^  which  are  constrained  to  have 
specific  values,  the  following  algorithm  solves  the  problem: 
maximize  f subject  to  the  constraints  on  the  variables  in  S . 

In  the  presentation,  we  do  not  distinguish  between  the  variables 
x^, ...,x^  and  the  corresponding  vertices  in  the  interaction  graph. 

Step  1.  If  n < 9 , solve  the  problem  by  exhaustively  trying  all 
possible  assignments  to  the  unconstrained  variables. 
Otherwise,  go  to  Step  2. 

Step  2.  Apply  Corollary  1 to  the  interaction  graph  G of  f . Let 
A,  B,  C be  the  resulting  vertex  partition.  Let  f^  b- 
the  restriction  of  f to  A pi  C and  let  f^  be  the 
restriction  of  f to  BljC  . For  each  possible  assigrunent 
of  values  to  the  variables  in  C-S  , perform  the  following 
steps : 

(a)  Maocimize  f^  with  the  given  vaulues  for  the  variables 
in  CUS  by  applying  the  method  recursively; 

(b)  Maximize  f^  with  the  given  values  for  the  variables 
in  CUS  by  applying  the  method  reciirsively ; 

(c)  Combine  the  solutions  to  (a)  and  (b)  to  obtain  a maximvim 
value  of  f with  the  given  values  for  the  variables 

in  C u S . 

Choose  the  assignment  of  values  to  variables  in  C-S  which 
maximizes  f and  return  the  ap|,ropriate  value  of  f as 


the  solution . 


The  correctness  this  algorithm  is  obvious.  If  n > 9 , the 

OfVtT ) 

algorithm  solves  at  most  2 _ subproblems  in  Step  2,  since  C 

is  of  0(Vn)  size.  Each  subproblem  contains  at  most 
2n/5  + 2 < 29n/50  variables.  Thus  if  t(n)  is  the  running 

time  of  the  algorithm,  we  have  t(n)  < 0(n  log  n)  + 2^^^^  • t(29n/50) 
if  n > 9 ^ ■t(n)  s 0(l)  if  n < 9 ' An  inductive  proof  shows 

that  t(n)  < 2^^"^  ^ . 
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h . Pebbling . 


The  following  one-person  game  arises  in  register  allocation 
problems  [21],  the  conversion  of  recursion  to  iteration  [l6],  and 
the  study  of  time-space  tradeoffs  [ it,10,l8].  Let  G = (V,E)  be 
a directed  acyclic  graph  with  maximum  in-degree  k . If  (v,  w")  is 
an  edge  of  G , v is  a predecessor  of  w and  w is  a 
successor  of  v . The  game  involves  placing  pebbles  on  the  vertices 
of  G according  to  certain  rules.  A given  step  of  the  game  consists 
of  either  placing  a pebble  on  an  empty  vertex  of  G (called  pebbling 
the  vertex)  or  removing  a pebble  from  a previously  pebbled  vertex. 

A vertex  may  be  pebbled  only  if  adl  its  predecessors  have  pebbles. 

The  object  of  the  game  is  to  successively  pebble  each  vertex  of  G 
(in  any  order)  subject  to  the  constraint  that  at  most  a given  number 
of  pebbles  are  ever  on  the  graph  simultajieously. 

It  is  easy  to  pebble  any  vertex  of  an  n-vertex  graph  in  n steps 
using  n pebbles.  We  are  interested  in  pebbling  methods  which  use 
fewer  than  n pebbles  but  possibly  many  more  than  n steps.  It  is 
known  that  any  vertex  of  an  n-vertex  graph  can  be  pebbled  with 
O(n/log  n)  pebbles  [10]  (where  the  constant  depends  upon  the  maximum 
in-degree),  and  that  in  general  no  better  bound  is  possible  [l8].  We 
shall  show  that  if  the  graph  is  planar,  only  0('/n)  pebbles  are 
necessary,  generalizing  a result  of  [l8].  An  example  of  Cook  [ 1 ] 
shows  that  no  better  bound  is  possible  for  planar  graphs. 

Theorem  3.  Any  n-vertex  planar  acyclic  directed  graph  with  maximm 
in-degree  k can  be  pebbled  using  0('/n  + k log^  n)  pebbles. 
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Proof.  Let  a = 2"j2  and  ^ = 2/5  . Let  G be  the  graph  to  be 

pebbled.  Use  the  following  recursive  pebbling  procedure.  If  n < n^  , 

idiere  n^  = (a/(l-p))  , pebble  all  vertices  of  G without  deleting 

pebbles.  If  n > n^  , find  a vertex  partition  A,  B,  C satisfying 

*/ 

Corollary  1.  Pesble  the  vertices  of  G in  topological  order.-' 

To  pebble  a vertex  v , delete  all  pebbles  except  those  on  C . For 
each  predecessor  u of  v , let  G(u)  be  the  subgraph  of  G induced 
by  the  set  of  vertices  with  pebble-free  paths  to  u . Apply  the  method 
recursively  to  each  G(u)  to  pebble  all  predecessors  of  v , leaving 
a pebble  on  each  such  predecessor.  Then  pebble  v . 

If  p(n)  is  the  maximum  number  of  pebbles  required  by  this  method 
on  any  n -vertex  graph,  then 

p(n)  = n if  n < n^  , 

p(n)  < a^Jn  + k + p(2n/5  + ot\fn  ) if  n > n^  . 

An  inductive  proof  shows  that  p(n)  is  0(  \/n  + k log^  n)  . □ 

It  is  also  possible  to  obtain  a substantial  reduction  in  pebbles 
while  preserving  a polynomial  boxand  on  the  number  of  pebbling  steps, 
as  the  following  theorem  shows. 

Theorem  U.  Any  n-vertex  planar  acyclic  directed  graph  with  maximum 
in-degree  k can  be  pebbled  using  0(n^^^  + k)  pebbles  in  0(kn^'^^)  time. 


That  is,  in  an  order  such  that  if  v is  a predecessor  of  w , 
V is  pebbled  before  w . 
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Proof.  Let  C be  a set  of  0(n^/^)  vertices  whose  removal  leaves 

*/  2/5 

G with  no  weakly  connected  component  containing  more  than  n ' 


vertices.  Such  a set  C exists  by  Theorem  2.  The  following  pebbling 
procedure  places  pebbles  permanently  on  the  vertices  of  C . pebble 
the  vertices  of  G in  topological  order.  To  pebble  a vertex  v , 


pebble  each  predecessor  u of  v and  then  pebble  v . To  pebble  a 
predecessor  u , delete  all  pebbles  from  G except  those  on  vertices 
in  C or  on  predecessors  of  v . Find  the  weakly  connected  component 
in  G minus  C containing  u . Pebble  all  vertices  in  this  component, 
in  topological  order. 

2/3 

The  total  number  of  pebbles  required  by  this  strategy  is  0(n  ' ) 

2/3 

to  pebble  vertices  in  C plus  n ' to  pebble  each  weakly  connected 
component  plus  k to  pebble  predecessors  of  the  vertex  v to  be 
pebbled.  The  total  number  of  pebbling  steps  is  at  most 
0(k.n.n^^^)  = 0(kn^/^)  . □ 


^ connected  component  of  a directed  graph  is  a connected 
component  of  the  undirected  graph  formed  by  ignoring  edge  directions. 
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5 . Lower  Bounds  on  Boolean  Circuit  Size. 


A Boolean  circuit  is  an  acyclic  directed  graph  such  that  each 
vertex  has  in-degree  zero  or  two,  the  predecessors  of  each  vertex  are 
ordered,  and  corresponding  to  each  vertex  v of  in-degree  two  is  a 
binary  Boolean  operation  b^  . With  each  vertex  of  the  circuit  we 
associate  a Boolean  function  ;Aiich  the  vertex  computes,  defined  as 
follows.  With  each  of  the  k vertices  v^  of  in-degree  zero  (inputs) 
we  associate  a variable  x.  and  an  identity  function  f (x.)  = x.  . 
With  each  vertex  w of  in-degree  two  having  predecessors  u , v we 
associate  the  function  f^  = b^(f^, f^)  . The  circuit  computes  the 
set  of  functions  associated  with  its  vertices  of  out-degree  zero 
(outputs). 

We  are  interested  in  obtaining  lower  bounds  on  the  size  (nimber 
of  vertices)  of  Boolean  circuits  which  compute  certain  common  and 
important  functions.  Using  Theorem  1 we  can  obtain  such  lower  bounds 
under  the  assumption  that  the  circuits  are  planar.  Any  circviit  can  be 
converted  into  a planar  circ\iit  by  the  following  steps.  First,  embed 
the  circuit  in  the  plane,  allowing  edges  to  cross  if  necessary.  Next, 
replace  each  pair  of  crossing  edges  by  the  crossover  circuit  illustrated 
in  Figure  1.  It  follows  that  any  lower  bound  on  the  size  of  planar 
circuits  is  also  a lower  bound  on  the  total  number  of  vertices  and 
edge  crossings  in  any  planar  representation  of  a non-planar  circuit. 

In  a technology  for  ^ich  the  total  niamber  of  vertices  and  edge 
crossings  is  a reasonable  meas\ire  of  cost,  o\ir  lower  bounds  imply  that 
it  may  be  expensive  to  realize  certain  commonly  used  functions  in 
hardware. 
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A superconcentrator  is  an  acyclic  directed  graph  with  m inputs 
and  m outputs  such  that  any  set  of  k inputs  and  any  set  of  k 
outputs  are  joined  by  k vert ex- disjoint  paths,  for  all  k in  the 
range  1 < k < m . 

Theorem  5 . Any  m- input,  m-output  planar  superconcentrator  contains 
2 

at  least  m /72  vertices. 


proof.  Let  G be  an  m- input,  m-output  planar  superconcentrator. 

Assign  to  each  input  and  output  of  G a cost  of  l/(2m)  , and  to  every 
other  vertex  a cost  of  zero.  Let  A,  B,  C be  a vertex  partition 
satisfying  Theorem  1 on  G (ignoring  edge  directions).  Suppose  C 
contains  p inputs  and  outputs.  Without  loss  of  generality,  suppose 
that  A is  no  more  costly  than  B , and  that  A contains  no  more 
outputs  than  inputs.  A contains  between  2m/5  - p and  m - p/2 
inputs  and  outputs.  Hence  A contains  at  least  m/5  - p/2  inputs 
and  at  most  m/2  - p/U  outputs.  B contains  at  least  m-p-  (m/2  - p/L) 
m/2  - 3p/^^  outputs.  Let  k = min{rm/5  - p/2l,rm/2  “ 3p/*+1  } • Since 
G is  a superconcentrator,  any  set  of  k inputs  in  A and  any  set  of 
k outputs  in  B are  joined  by  k vert ex- disjoint  paths.  Each  such 
path  must  contain  a vertex  in  C vdiich  is  neither  an  input  nor  an 
output.  Thus  - p > min(m/3  - p/2,  m/2  - 5p/*+}  > m/3  " P f 

and  n > m /72  . □ 


The  property  of  being  a superconcentrator  is  a little  too  strong 
to  be  useful  in  deriving  lower  bo\inds  on  the  complexity  of  interesting 
functions.  However,  there  are  weaker  properties  which  still  require 

p 

fi(m  ) vertices.  Let  G = (V,  E)  be  an  acyclic  directed  graph  with  m 


numbered  inputs  v-,v  , and  m numbered  outputs  w,,w„,  ...,w  . 

X c m X 2 IR 

G is  said  to  have  the  shifting  property  if,  for  any  k in  the  range 
1 < k < m , ajiy  I in  the  range  0 < I < m-k  , and  any  subset  of  k 


soTirces  {v.  } such  that  i, , i„,  . . i,  < m-/ 

il  d K- 

vertex-dis Joint  paths  joining  the  set  of  inputs  [v.  , 

^1 

the  set  of  outputs  +i'*"’''^i  +l^  ' 


there  are  k 
.,v^  } with 


Theorem  6.  Let  G be  a planar  acyclic  directed  graph  with  the 

O 

shifting  property.  Then  G contains  at  least  L®/2J  /l62  vertices. 

Proof.  Suppose  that  G contains  n vertices.  Assign  a cost  of  l/m 
to  each  of  the  first  [_m/2 ^ inputs  and  to  each  of  the  last  Lm/2J 
outputs  of  G , and  a cost  of  zero  to  every  other  vertex  of  G . Call 
the  first  inputs  and  the  last  outputs  of  G costly. 

Let  A,B,C  be  a vertex  partition  satisfying  Theorem  1 on  G 
(ignoring  edge  directions). 

Without  loss  of  generality,  suppose  that  A is  no  more  costly 
than  B , and  that  A contains  no  more  costly  outputs  than  costly 
inputs.  Let  A'  be  the  set  of  costly  inputs  in  A , B'  the  set  of 
costly  outputs  in  B , p the  number  of  costly  inputs  and  outputs 
in  C , and  q the  number  of  costly  inputs  and  outputs  in  A . Then 
2LV2j/3  - P < q < LV2J  - p/2  . Hence  jA’ | > q/2  > Lm/2j/3  - p/2  . 
Also 
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lA-  1-  IB-  1 > lA'  l-(Lm/2j  - P - (q-  lA'  |)) 
> q/2  • (Lm/2j  - p - q/2) 


> (Lni/2j/3  - p/2)(Lm/2J  - p - + p/2) 

= (LV2J/5  - p/2)(2Lm/2j/3  - p/2) 

> 2Lm/2j^/q  - pLm/2j/2  . 


For  V.  e A'  , w.  e B’  , and  l in  the  range  1 < ■*  < L“/2j  > 

call  , e a match  if  j-i  = f . For  every  e A'  and 

w.  cB'  there  is  exactly  one  value  of  l which  produces  a match; 

J 

hence  the  total  number  of  matches  for  all  possible  v.  , w.  , / is 

A J 

p 

Ia'1'IB'I  > 2Lm/2j  /q  - pj_m/2j/2  . Since  there  are  only  Lm/2j 

values  of  I , some  value  of  I produces  at  least  2Lm/2j/q  - p/2 

matches.  Thus,  for  k = 2Lm/2j/q  - p/2  , there  is  some  value  of  I 

and  some  set  of  k inputs  A"  = (v.  ,v.  ,...,v.  } c A'  such  that 

^1  2 k 

B"  = (w.  } c B'  . Since  G has  the  shifting 

property,  there  must  be  k vertex-disjoint  paths  between  A"  and  B"  . 
But  each  such  path  must  contain  a vertex  of  C which  is  neither  an 
input  nor  an  output.  Hence  2^J2’Jn  -p  > 2|_m/2j/q-  p/2  , and 
n > Lm/2j^/l62  . □ 


A shifting  circuit  is  a Boolean  circuit  with  m primary  inputs 
^l’^2’*'*’^m  ’ or  more  aiixiliary  inputs,  and  m outputs  z^,  z^,  ...,z^  , 

such  that,  for  any  k in  the  range  0 < k < m , there  is  some  assignment 
of  the  constants  0,1  to  the  auxiliary  inputs  so  that  output 
computes  the  identity  function  x^  , for  ^ • Th®  Boolean 


L 


It. 


I 


convolution  oi  two  Boolean  vectors  (x, ) and  (v, , v-i . . . . v ) 
12  m ''m 

is  the  vector  z^, . . . , given  by  z^  = Xj  . 

i+J  = k 


• 2 
Corollary  2.  Any  planar  shifting  circuit  has  at  least  L’^i/2J  /l62 

vertices. 

Proof.  Any  shifting  circuit  has  the  shifting  property. 

See  [25,2k].  □ 

Corollary  5«  Any  planar  circuit  for  computing  Boolean  convolution  has 
at  least  L^/^J  vertices. 


Proof.  A circuit  for  computing  Boolean  convolution  is  a shifting 

circuit  if  we  regard  x^,...,x  as  the  primary  inputs  and  z„, ...,z 

i-  K1  2 ITt+'l 

as  the  outputs,  □ 


Corollary  k.  Any  planar  circuit  for  computing  the  product  of  two  m 

2 

bit  binary  integers  has  at  least  Liii/2J  /l62  vertices. 


Proof.  A circuit  for  m\altiplying  two  m-bit  binary  integers  is  a 
shifting  circiait.  □ 


The  last  result  of  this  section  is  an  Q(m  ) lower  bound  on  the 
size  of  any  planar  circuit  for  miiltipiying  two  mxm  Boolean  matrices. 
We  shall  assvune  that  the  inputs  are  x.  . , y . . for  1 < i,  j < m and 
the  outputs  are  z^j  for  1 < i,J  < m . The  circuit  canputes 
Z = X*Y  , where  Z = (z^^)  , X = (x^^)  , and  Y = (y^j)  . We  use 
the  following  property  of  circuits  for  multiplying  Boolesin  matrices, 
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called  the  matrix  concentration  property  [2^,2h].  For  any  k in  the 

range  1 < k < n , any  set  {x.  . | 1 < r < k}  of  k inputs  from  X , 

ir^r 

and  any  permutation  o of  the  integers  one  throvigh  n , there  exist 

k vertex-disjoint  paths  from  '(x.  . 1 1 < < k}  to  [z.  / . v | 1 < r < k] 

^r*^r  2 ''^r-  ~ 

Similarly,  for  any  k in  the  range  1 < k < n , any  set 

t^i  .1  1 1 < r < k}  of  k inputs  from  Y , and  any  permutation  cr 

of  one  throu^  n , there  exist  k vertex-disjoint  paths  from 

[y±  , 1 1 r < k]  to  [z^,  V 1 1 < r < k]  . 


Theorem  7»  Any  planar  circuit  G for  multiplying  two  m x m Boolean 
matrices  contains  at  least  cm  vertices,  for  some  positive  constant  c 


Proof.  This  proof  is  some^at  involved,  and  we  make  no  attempt  to 

maximize  the  constant  factor.  Suppose  G contains  n vertices,  and 

2 

that  m is  even.  Assign  a cost  of  l/(Um  ) to  each  input  x^^  and 

2 

each  input  y, . , a cost  of  l/(2m  ) to  each  output  z.  . , and  a cost 
ij  ^ J 

of  zero  to  every  other  vertex.  There  is  a partition  A , B , C of  the 
vertices  of  G such  that  neither  A nor  B has  toteuL  cost  exceeding 
1/2  , no  edge  joins  a vertex  in  A with  a vertex  in  B , and  C 
contains  no  more  than  2’j2'jn  / (1  - •'/2/3 ) = c^4n  vertices.  This 
is  a corollary  of  Theorem  1;  see  [lU  ] . Without  loss  of  genersQ-ity, 
suppose  that  B contains  no  fewer  outputs  than  A , and  that  A 
contains  no  fewer  inputs  than  inputs  y^^  . Then  B contains 

at  least  (m^  - c^^Vn  )/2  outputs,  which  contribute  at  least 
l/U  - Cj^'/rT/ (Um^)  to  the  cost  of  B . Thus  inputs  contribute  at  most 
l/h  - Cj^Vn/(Um^)  to  the  cost  of  B , and  B contains  at  most 


I 


F 


n*"  *■  c^-/n  inputs.  A contains  at  least  2;n'  - (in'  t c^  fn  ) - c^Vn  = 
n'"  - 2c^'^  inputs,  of  which  at  least  rti~ /2  - c^  \fn  are  inputs  . 

.Tne  of  the  fol Lowing  cases  must  hold. 


Case  1.  A contains  at  least  3in  /5  inputs  x,.  ^ . Let  p he  the 

^ J 

number  of  columns  of  X ’A/hich  contain  at  least  hm/'J  elements  of  A . 

o 

Then  pm  (m-p)(Um/7)  > 3m^/5  , and  p > m/l5  . Let  q be  the  number 

of  columns  of  Z which  contain  at  least  l»m/9  elements  of  B . Then 

qm+  (m-q)(l+m/q)  > m"/2  - c^-Jn /2  , and  q > m/lO  - 9c^*/n /(lOm)  . 

Let  k = min{m/l5  , m/lO  - 9c^'iii/(10m)]  . Choose  any  k columns 

of  X , each  of  which  contains  at  least  l+m/7  elements  of  A . Match 

each  such  col\imn  of  X with  a column  of  Z which  contains  at  least 

Um/9  elements  of  B . For  each  pair  of  matched  columns  x„.  , z„  . , 

' *1  *0 

select  a set  of  W/7  + W/9-m  = m/63  rows  t such  that  x^^  is 

in  A and  z . is  in  B . Such  a selection  gives  a set  of  km/63 
^ J 

elements  in  XOA  and  a set  of  km/63  elements  in  ZflB  vAiich  must 
be  joined  by  km/63  vertex-disjoint  paths,  since  G has  the  matrix 
concentration  property.  Each  such  path  must  contain  a vertex  of  C . 

Thus  km/63  < c^\/n  , which  means  either  m^/(15*63)  < c^-Jn  (i.e., 
(m^/(l5*65c2_))^  < n ) or  m/63(m/lO  - 9c^'\/n  / (10m) ) < c^\In 
(i.e.,  (m^/(9*69c^))^  < n ). 

O 

Case  2.  A contains  fewer  than  3m  /?  inputs  x^^  . Then  A contains 

at  least  2m  - 2o  vn  inputs  y.  . . Let  S be  the  set  of  m/2 

1 J 

columns  of  Z which  contain  the  most  elements  in  B . 


19 


Subcase  2a. 


p 

S contains  at  least  /lO  elements  in  B . Let  p 


be  the  number  of  coliamns  of  X which  contain  at  least  km/9  elements  j 

of  A . Then  pm+U(m-p)m/9  > m^/2  - c^VtT  , and  p > m/lO  - 9Cj^Vn/(5m)  . 

Let  q be  the  number  of  columns  of  Z vdiich  contain  at  least  km/j 

p 

elements  of  B . Then  qm+  U(m/2  - q)m/7  > 3m“/l0  , and  q > m/jO  . 

A proof  similar  to  that  in  Case  1 shows  that  n > cm  for  some  positive 
constant  c . 

Subcase  2b.  S contains  fewer  than  / 10  elements  in  B . Then  the 

m/2  col\jmns  of  Z not  in  S contain  at  least  m^ /5  - /2  elements  * 

in  B . Let  q be  the  number  of  columns  of  Z not  in  S vdiich  contain 
at  least  m/lO  elements  in  B . Then  (5n+  (m/2  - q)  (m/lO)  ^ ro^/5  - c^Vn'/2  , 

and  q > m/6  - 5c^'\/rr/ (9m)  . If  0 > q _>  m/6  - 5Cj^Vn/(9ni)  > then  / 

2 2 

(3m  /(lOc^))  > n . Hence  assume  q > 0 . Then  all  columns  in  S . 

must  contain  at  least  m/lO  elements  in  B , and  2m/3  - 5c^V^/(9m) 
columns  of  Z must  contain  at  least  m/lO  elements  in  B . 

Let  p be  the  number  of  col\mins  of  Y which  contain  at  least  m/25 
elements  of  A . Then  pm+  (m-p)(m/25)  2m^/5  - 2c^Vn  , and 

p > 3m/8  - 25c^'\/n /(12m)  . 

For  any  input  y. . e A and  integer  t in  the  range  -n+1  < t < n-1  , 

call  y..  I a match  if  ^ e B . By  the  previous  con5)utations, 

ij  

there  are  at  least  2m/3  - 5Cj^'/n  / (9®)  + 5®/8  - 23c^'/n  / (12m)  - m s= 
m/25  - 95cj^'\/ir/(36m)  = m/25  - c^VrT/m  colianns  j such  that  y^j 

contains  m/25  elements  of  A and  contains  m/lO  elements  j 

of  B . Each  such  column  produces  m /250  matches;  thus  the  toteLL 
nxrnber  of  matches  is  at  least  m^ /6250  - mc^Vrr/250  . Since  there  are 
only  2m-l  values  of  I , some  value  of  i produces  at  least 


k = in^/12, 500  - /500  matches.  Since  G has  the  matrix 

concentration  property,  this  set  of  matches  corresponds  to  a set 
of  k elements  in  YQA  and  a set  of  k elments  in  ZOB  which 
must  he  joined  by  k vertex-disjoint  paths.  Each  such  path  must 
contain  a vertex  in  C . Thus  k < c^Vn  , which  means 
mV(12,500(c^  + C2/500))^  < n . 

4 

In  all  cases  n > cm  for  some  positive  constant  c . Choosing 
the  minimian  c over  all  cases  gives  the  theorem  for  even  m . The 
theorem  for  odd  m follows  inmediately.  □ 

The  bo\inds  in  Theorems  5-7  and  Corollaries  2 - U are  tight  to 
within  a constant  factor.  We  leave  the  proof  of  this  fact  as  an 
exercise. 


I 


21 


6.  Embedding  of  Data  Structures. 

Let  = (V^, and  = {^21^2^  136  \indirected  graphs.  An 

embedding  of  G^^  in  Gg  is  a one-to-one  map  0:  -*  Vg  • The  vorst-case 

proximity  of  the  embedding  is  max{d2(^(v),)6(w) ) | [v^wjeE^}  , where 
d2(x,y)  denotes  the  distance  between  x and  y in  G^  • The  average 
proxi mi ty  of  the  embedding  is  E {d2(0(v)^ji(w))  1 {v,w}eE^}  . 

These  notions  arise  in  the  following  context.  Suppose  we  wish  to 
represent  some  kind  of  data  structure  by  another  kind  of  data  structure, 
in  such  a way  that  if  two  records  are  logically  adjacent  in  the  first 
data  struct\ire,  their  representations  are  close  together  in  the  second. 

We  can  model  the  data  structiires  by  \mdirected  graphs,  with  vertices 
denoting  records  and  edges  denoting  logical  adjacencies.  The  representation 
problem  is  then  a graph  embedding  problem  in  which  we  wish  to  minimize 
worst-case  or  average  proximity.  See  [5,13»19]  Tor  research  in  this  area. 

Theorem  8.  Any  planar  graph  with  maximum  degree  k can  be  embedded  in 
a bineiry  tree  so  that  the  average  proximity  is  a constant  depending  only 
upon  k . 


Proof.  Let  G be  sui  n- vert ex  planar  graph.  Embed  G in  a binary 

tree  T by  using  the  following  reciirsive  procedure.  If  G has  one 
vertex  v , let  T be  the  tree  of  one  vertex,  the  image  of  v . 
Otherwise,  apply  Corollary  1 to  find  a partition  A , B , C of  the 
vertices  of  G . Let  v be  any  vertex  in  C (if  C is  enqjty,  let  v 
be  any  vertex).  Qnbed  the  subgraph  of  G induced  by  AuC-{v}  in  a 
binary  tree  T^^  by  applying  the  method  recursively.  Embed  the  subgraph 
of  G induced  by  B in  a binary  tree  Tg  by  applying  the  method 
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Let  h(n)  be  the  maxim-um  depth  of  a tree  T of  n vertices 
produced  by  this  algorithm.  Then 


h(n)  < 9 if  n < 9 , j 

h(n)  < h(2n/5  + aViVn  - 1)  < h(29ti/50)  if  n > 9 . 

It  follows  that  h(n)  is  O(log  n)  . 

Let  G = (V,  E)  be  an  n-vertex  graph  to  vdiich  the  algorithm  is  I ■ 

applied,  let  be  the  subgraph  of  G induced  by  AIJC  , and  let 

G^  be  the  subgraph  induced  by  B . If  s(G)  = {d.2(0(v),0(w) ) [ (v,w)  e E]  , 

then  s(G)  = 0 if  n = 1 , and  s(G)  < s(G^)  + s(G2)  + k |C  | h(n)  if  ^ 

n > 1 . This  follows  from  the  fact  that  any  edge  of  G not  in  G^  or  G^ 
must  be  incident  to  a vertex  of  C . 

If  s(n)  is  the  maximum  value  of  s(g)  for  any  n-vertex  graph  G , 

then 

s(l)  = 0 ; 

s(n)  < max{s(i)  + s(n-i)  + ckVn  log  n[  n/5  - < i < 2n/5  + 2'j2\fn  } 

if  n > 1 , for  some  positive  constant  c . ■ 

An  inductive  proof  shows  that  s(n)  is  0(kn)  . 

If  G is  a connected  n-vertex  graph  embedded  by  the  algorithm,  then 
G contains  at  least  n-1  edges,  and  the  average  proximity  is  0(k)  . 

If  G is  not  connected,  embedding  each  connected  component  separately 

> 

and  combining  the  resulting  trees  arbitrarily  achieves  an  0(k)  average 
proximity.  □ 
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i i 

< It  iij  natural  to  ask  whether  any  graph  of  bounded  degree  can  be 

I j 

I ; embedded  in  a binary  tree  with  0(l)  average  proximity.  (Graphs  of 

unbounded  degree  cannot  be  so  embedded;  the  star  of  Fig\u:e  2 requires 
•il(log  n)  proximity.)  Such  is  not  the  case,  and  in  fact  the  proj'erty 
of  being  embeddable  in  a binary  tree  with  0(l)  average  proximity  is 
closely  related  to  the  property  of  having  a good  separator. 

To  make  this  statement  more  precise,  let  S be  a class  of  graphs. 

The  class  S has  an  f(n)  -separator  theorem  if  there  exist  constants 
f a 1 , p 0 such  that  the  vertices  of  any  n-vertex  graph  in  S can 

I be  } artitioned  into  three  sets  A,  B,  C such  that  |a|,  |b|  < an  , 

! I^i  ^ no  vertex  in  A is  adjacent  to  any  vertex  in  B . 

Let  S be  any  class  of  graphs  of  bounded  degree  closed  under  the 
subgraph  relation  (i.e.,  if  Gg  e S and  G^^  is  a subgraph  of  G^  , 

I then  G^ c S ).  Suppose  S satisfies  an  ng(n)/(log  n)  separator 

’ theorem  for  some  non-decreasing  function  g(n)  . Using  the  idea  in 

' the  proof  of  Theorem  8,  it  is  not  hard  to  show  that  any  graph  in  S 

can  be  embedded  in  a binary  tree  with  0(g(n))  average  proximity. 
Conversely,  suppose  any  graph  in  a class  S can  be  embedded  in  a binary 
tree  with  0(g(n))  average  proximity.  Then  S satisfies  an  ng(n)/log  n 
separator  theorem.  In  particular,  if  S satisfies  no  o(n)  -separator 
I theorem,  then  embedding  the  graphs  of  S in  binary  trees  requires 

Q(log  n)  average  proximity.  ErdOs,  Graham,  and  Szemeredi  [ 7 ] have 
shown  the  existence  of  a class  of  graphs  of  bounded  degree  having  no 
o(n)  -separator  theorem. 


i 

I 
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7.  The  Post  Office  Problem. 


In  [U],  laiuth  mentions  the  following  problem:  given  n points 
(post  offices)  in  the  plane;  determine,  for  any  new  point  (house), 
which  post  office  it  is  nearest.  Any  preprocessing  of  the  post  offices 

is  allowed  before  the  houses  ar_  processed.  Shamos  [22]  gives  an 

2 2 
O(log  n)  -time,  0(n  ) -space  algorithm  and  an  O((log  n)  ) -time, 

0(n  log  n)  -space  algorithm.  See  also  [6],  Using  Theorem  2 we  can 

give  a soltition  which  requires  O(log  n)  time  and  0(n)  space,  both 

minimum  if  only  binary  decisions  are  allowed, 

A polygon  is  a connected,  open  planar  region  bounded  by  a finite 
set  of  line  segments.  (For  convenience,  we  allow  the  point  at  infinity 
to  be  an  endpoint  of  a line  segment;  thus  a line  is  a line  segment.) 

A polygon  partition  of  the  plane  is  a partition  of  the  plane  into 
polygons  and  bounding  line  segments.  A tri angulation  of  the  plane  is 
a polygon  partition,  all  of  \diose  polygons  are  bomded  by  three  line 
segments.  A tri angulation  of  a polygon  partition  is  a refinement  of 
the  partition  into  a triangulation.  Two  polygons  in  a polygon  partition 
are  adjacent  if  their  boTondaries  share  a line  segment.  A set  of  polygons 
is  connected  if  any  two  polygons  in  the  set  are  joined  by  a sequence  of 
adjacent  polygons. 

We  shall  solve  the  following  triangle  problem:  given  an  n-triangle 
triangulation  and  a point,  determine  vAiich  triangle  or  line  segment  of 
the  triangulation  contains  the  point.  The  post  office  problem  can  be 
reformulated  as  a triangle  problem;  the  set  of  points  closest  to  each 
post  office  forms  a polygon  [22],  We  shall  make  use  of  the  following 
lemma,  which  we  do  not  prove. 
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Lemma  1.  Any  n-polygon  partition  has  a refinement  ^ose  total  number 


of  triangles  is  bovuided  by  n plus  the  number  of  line  segments  bounding 
non-triangles  plus  a constant  (a  line  segment  bounding  two  non-triangles 
counts  twice  in  this  bound). 

We  shall  b-uild  up  a sequence  of  more  and  more  con^jlicated  (but 
more  and  more  efficient)  algorithms,  the  last  of  which  is  the  desired  one. 

Theorem  9«  Given  an  O(log  n)  -time,  0(n^^^)  -space  algorithm  for  the 
triangle  problem  with  e > 0 , one  can  construct  an  O(log  n)  -time, 

-space  algorithm. 

Proof.  Let  T be  a triangulation  and  v be  a vertex  for  viiich  the 
triangle  problem  is  to  be  solved.  By  Theorem  2 there  is  a set  of  0(n^/^) 
triangles  Cq  whose  removal  from  T leaves  no  connected  set  of  more  than 
0 (n^/^ ) t ri angle  s . 

Merge  pedrs  of  adjacent  triangles  which  are  not  in  Cq  to  form  a 

2/3 

polygon  partition  Pq  . Pq  contains  at  most  0(n  '■^)  line  segments, 
since  each  such  line  segment  must  be  a bounding  segment  of  a triangle 
in  T . Find  a triangulation  Tq  of  Pq  with  0(n^/^)  triangles, 
vrtiich  exists  by  Lemma  1.  Using  the  given  algorithm,  determine  >Aiich 
triangle  or  line  segment  of  Tq  contains  v . 

If  V is  in  some  triangle  of  Cq  , the  problem  is  solved.  Otherwise, 

V is  known  to  be  in  some  connected  set  C^^  of  triangles  in  T minus  Cq  . 
Merge  pairs  of  adjacent  triangles  which  are  not  in  C^  to  form  a polygon 
partition  P^  . Since  P^  contains  at  most  0(n  ' ) line  segments, 
there  is  a triangulation  T^^  of  P^^  with  triangles.  Using 

the  given  CLlgorlthm,  determine  which  triangle  or  line  segment  of  T^ 
contains  v . This  solves  the  problem. 
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The  sets  , polygon  partitions  , and  triangulations 

are  all  precomputed.  Thus  the  time  required  by  the  algorithm  is 

O(log  n^/^)  to  discover  \diich  triangle  of  Tq  contains  v , plus 

O(log  n^/^)  to  discover  \diich  triangle  of  T^  contains  v . The 

total  time  is  thus  O(log  n)  . The  total  space  is 

T.  OClT^l^®)  < . □ 

i 

Corollary  For  any  e > 0 there  is  an  O(log  n)  -time,  0(n^^^) 

-space  algorithm  for  the  triangle  problem. 

Proof.  Immediate  from  Theorem  9>  using  the  O(log  n)  -time, 

O 

0(n  ) -space  algorithm  of  [22]  as  a starting  point.  □ 

Theorem  10.  There  is  an  0(log  n)  -time,  0(n)  -space  solution  to  the 
triangle  problem. 

Proof.  Let  T be  a triangulation  and  v a vertex  for  which  the  triangle 
problem  is  to  be  solved.  If  T contains  no  more  than  n^  triangles, 
where  n^  is  a sufficiently  large  constant,  determine  -vAiich  triangle 
contains  v by  testing  v against  each  line  segment  boiinding  a triangle 
of  T . Otherwise,  let  C be  a set  of  0(n^/^)  triangles  whose  removal 
from  T leaves  no  connected  set  of  more  than  0(n^'^^)  triangles.  Group 
the  connected  sets  of  triangles  in  T minus  Cq  into  sets  , each 

V5 

containing  within  a constant  factor  of  n ' triangles. 

Merge  pairs  of  adjacent  triangles  which  are  not  in  Cq  to  form 

3 /5 

a polygon  partition  Pq  . Pq  contains  at  most  0(n^'  ) line  segments. 
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Find  a triangulation  Tq  of  with  0(n^/^)  triangles.  Using  an 

O(log  n)  -time,  0(n'^'^^)  -space  algorithm,  determine  vdiich  triangle 
of  Tq  contains  v . 

If  V is  some  triangle  of  Cq  , the  problem  is  solved.  Otherwise 
V is  known  to  be  in  some  set  . Merge  pairs  of  adjacent  triangles 
which  are  not  in  to  form  a polygon  partition  . Each  line 

segment  bounding  a non- triangular  polygon  of  P^^  must  bound  a triangle 
of  Cq  . Thus  there  is  a tri angulation  of  P^  containing 

|C^|  + 0(n^/^)  triangles.  Apply  the  algorithm  recursively  to  discover 
which  triangle  of  contains  v . This  solves  the  problem. 

The  sets  , polygon  partitions  P^^  , and  triangulations  T^ 
are  all  precomputed.  If  t(n)  is  the  worst-case  time  required  by  the 
algorithm  on  an  n-triangle  triangulation,  then 
t(n)  = 0(1)  if  n < nQ  , 

t(n)  = t(0(n^'^^ ) ) + 0(log  n)  otherwise. 

An  inductive  proof  shows  that  t(n)  is  0(log  n)  if  Uq  is  chosen 
sufficiently  large. 

If  s(n)  is  the  worst-case  storage  space  required  by  the  algorithm 
on  an  n-triangle  tri angulation,  then 
s(t)  = 0(1)  if  n < Uq  , 

s(n)  < 0(n^''^^^)  + max{Z/  s(n^  + 0(n^/^))  | Z!  <n  and 

< n^^  < c^n^/^} 

for  some  positive  constants  c^  and  c^  . 

An  inductive  proof  shows  that  s(n)  is  0(n)  . □ 
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The  preprocessing  time  required  by  the  algorithm  of  Theorem  10 
is  0(n  log  n)  . See  [22],  We  do  not  advocate  this  algorithm  as  a 


practical  one,  but  its  existence  suggests  that  there  may  be  a practical 
algorithm  with  an  0(log  n)  time  bound  and  an  0(n)  space  bound. 
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8.  other  Applications. 


As  illustrated  in  this  paper.  Theorem  1 and  its  corollaries  have 
many  interesting  applications,  and  the  paper  does  not  exhaust  them. 

We  have  obtained  two  additional  resiilts  \rtiich  require  fuller  discussion 
than  is  possible  here.  One  is  the  application  of  Theorem  1 to  Gaussian 
elimination.  George  [8]  has  proposed  an  0(n  log  n)  -space,  0(n^'^^)  -time 
method  of  carrying  out  Gaussian  elimination  on  a system  of  equations  whose 
sparsity  structure  corresponds  to  a ^/n  x Vn  square  grid.  We  can 
generalize  his  method  so  that  it  applies  to  any  system  of  equations 
whose  sparsity  structure  corresponds  to  a planar  or  almost-planar  graph. 
Such  systems  arise  in  the  solution  of  two-dimensional  finite -element 
problems  [15].  We  shall  discuss  this  application  in  a subsequent  paper; 
we  hope  that  it  will  prove  of  practical,  as  well  as  theoretical,  value. 

Another  application  involves  the  power  of  non-determinism  in  one-tape 
Turing  machines.  We  can  prove  that  any  non-deterministic  t(n)  -time- 
bounded  one-tape  Turing  machine  can  be  simulated  by  a t(n)^  alternating 
one-tape  Turing  machine  with  a constant  number  of  alternations,  where 
> < 1 is  a suitable  constant  and  t(n)  satisfies  certain  reasonable 
restrictions.  Alternation  generalizes  the  concept  of  non- determinism 
and  is  discussed  in  [5,12].  Our  result  strengthens  Paterson's  space- 
efficient  simulation  of  one-tape  Turing  machines  [I7]. 
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